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Abstract 

The effective macroscopic response of nonlinear elastomeric inhomogeneous 
materials is of great interest in many applications including nonlinear com- 
posite materials and soft biological tissues. The interest of the present work 
is associated with a microsphere composite material, which is modelled as a 
matrix-inclusion composite. The matrix phase is a homogeneous isotropic 
nonlinear rubber-like material and the inclusion phase is more complex, con- 
sisting of a distribution of sizes of stiff thin spherical shells filled with gas. 
Experimentally, such materials have been shown to undergo complex defor- 
mation under cyclic loading. Here, we consider microspheres embedded in an 
unbounded host material and assume that a hydrostatic pressure is applied 
in the 'far-field'. Taking into account a variety of effects including buck- 
ling of the spherical shells, large deformation of the host phase and evolving 
microstructure, we derive a model predicting the pressure-relative volume 
change load curves. Nonlinear constitutive behaviour of the matrix medium 
is accounted for by employing neo-Hookean and Mooney-Rivlin incompress- 
ible models. Moreover a nearly-incompressible solution is derived via asymp- 
totic analysis for a spherical cavity embedded in un unbounded isotropic 
homogeneous hyperelastic medium loaded hydrostatically. The load-curve 
predictions reveal a strong dependence on the microstructure of the compos- 
ite, including distribution of microspheres, the stiffness of the shells, and on 
the initial volume fraction of the inclusions, whereas there is only a modest 
dependence on the characteristic properties of the nonlinear elastic model 
used for the rubber host. 

Keywords: microsphere, composite, pressure- volume curve, buckling, 
nonlinear elasticity, rubber, Mooney-Rivlin 



1. Introduction 



Microsphere composites are used in a multitude of industrial applica- 
tions. Good examples are ultra-low density fillers in engineering materials 
such as composites, coatings, sealants, explosives, automotive components, 
paint and crack filler s and elastomers and as blowing agents in printing inks 
flAsh and Ashl . l2007l ). The microsphere is a spherical particle, a few mi- 
crons in diameter, with a thermoplastic shell, and shell to diameter ratio 
typically of the order of 0.01. The use of microspheres in materials brings 
forth numerous benefits which include reducing density, improving stability, 
increasing impact strength, providing a smoother surface finish, increasing 
thermal insulation, increasing compressibility and often reducing costs. A 
scanning electron microscopy image of the microstructure of a hollow glass 
microspher e compos i te (w ith very high volume fraction) is shown in Fig. [1] 
taken from iLi et al.l ( 1201 ll ). 




(a) 



(b) 



Figure 1: Scanning electron microscopy image of (a) a Hollow Glass Microsphere composite 
material and (b) its surface/shell structure. In this situation (with applications to thermal 
conductivity) the compos i te is d ensely filled with microspheres. Reproduced (with kind 
permission) from lLi et al.l ( 201l[ ). 



The application of interest in this paper is that of acoustics, using mi- 
crosphere composites as a means of reducing sound reflection. The com- 
posite consists of an elastomeric matrix phase, inside which are located a 
large number of ran domly distributed Expancel microspheres, see Fig 
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Shorter et al.l ( 120081 ). Of specific interest is how sound reflection can be 
affected by a macroscopic hydrostatic pressure applied to the microsphere 
composite. These materials have been found to be useful in such conditions 
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because the presence of the reinforcing shell delays the onset of the cavity 
collapse and the consequent degradation in the acoustic performance of the 
composite. In order to understand exactly how the acoustic characteristics of 
the material are affected by applied pressure, it is necessary to develop mod- 
els that describe how the composite deforms mechanically under this applied 
loading. Experimentally it is known that the constitutive pressure-relative 
volume change curve is nonlinear (and hysteretic during unloading) but the 
dominant physical mechanisms contributing to this nonlinearity are still not 



fully understood. In Fig. 2(b) we illustrate the constitutive behaviour of the 
material with some experimentally determined load curves associated with 
the composite for i ncreasing volume fractions of the microsphere material 
Shorter et all fcoosh . 
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Figure 2: In (a) we show an image of a Silicone Room Temperature Vulcanizing (RTV) 
microsphere elastomer filled with 5% volume of Expancel microspheres. In (b) we show 
an experimentally determined stress-strain curve associated with this composite under 
uniaxial tension. The solid curve is unfilled and others refer to incre asing volume fraction s 
of the microsphere phase. Reproduced (with kind permission) from lShorter et aL ( 2008 ). 



A wide range of work on the modelling of microsphere filled composites 
has appeared in the literature. For the prediction of their acoustic proper- 
ties without applied pressure, a ii umbe r of models have b een proposed, see 
e.g. iGaunaurd and Uberalll ( 1l982l ) and iBaird et al.l (119991 ) in the elastic and 
viscoelastic case respectively. Such predictions typically rely upon the use of 
classical multiple scattering models, the most commonly used being those 
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of IWaterman and Truelll (119611 ) : iKuster and Toksozl (1197411: iBose and Ma] 



(Il974h : iGaunaurd and Uberalll (Il982h : lAnson and Chiverg (jl989l ). We note 
th at a useful compari s on of experiments and various theories can be found 



m 



Anson and ChiversI (Il993[) . 



Gaunaurd et all (119841 ) considered the pressure dependence of dynamic 



moduh, albeit in a simplified case of porous solids, i.e. in the absence of the 
shell phase and thus the effect of the pressure is to reduce the pore size, 
the main interest lying in the dynamic material properties. Here, although 
we are certainly interested in the dynamic response, we wish first to un- 
derstand the origins of nonlinearity in the pressure-relative volume change 
loading curve associated with the composite. Significant work has be e n car- 
ried out on the deformation of porous materials, see e.g. [Mackenzie! (119501 ) 
for an early model for associa ted effective linear elastic properties, and for 
rubbe r foar n materials, s ee e.g . iGent and ThomasI (jl959[ ). iGibson and Ashby 
( 119821 ) and iLakes et al.l ( 1l993l ) , where the principal mechanisms of deforma- 
tion are well understood. However, the composite under consideration here 
has a more complex microstructure due principally to the presence of mi- 
crospheres and the lack of understanding as to how they behave within the 
elastomeric substrate under apphed pressure. Microspheres are also present 
in the context of syntactic foams (a material comprising a polymeric matrix 
filled with microspheres). Such materials open up the possibility of low den- 
sity materials with high tolerance to damage. Much of the recent modelling 
work in this area however has focused on the effective linear elastic properties 



used, see e.g. 


Bardella and Genua ( 


2001) 


Porfiri and Gupta 


(2009h: 


Taeliavia et al. 



(120091). 



Few models deal with the nonlinea r response of a microsphere composite 
under loading. iKerr and BairdI (j2002[ ) proposed an elasto-plastic model for 
the load curve and subsequent prediction of dynamic material properties. 
One criticism of this model would be that plasticity yields permanent de- 
formation. However, it is well acknowledged that although the load-unload 
curve is hysteretic, when all load has been removed the rn aterial (eventually) 
returns to its original configuration (IBrazier-Smithl . 120101 ). Therefore it does 
not appear that plastic deformat i on is the cause of nonlinearity. In a re- 
lated application iPanigrahi et al.l (|2008[ ) considered the acoustic response of 
inhomogeneous media under applied pressure, although it ap pears that the 



microstructure is rather different from that considered here. In Shorter et al 



( 120081 ) ■ [Shorter et al.l ( l2010l ) the problem of the buckling of a single, isolated 
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spherical shell was considered (as a result of axial compression, rather than 
hydrostatic pressure) using Finite Element Analysis and results were subse- 
quently compared with experiments involving a table tennis ball embedded 
inside a transparent elastomer. Comparisons were made between perfectly 
bonded and unbonded spherical shells and subsequent buckling response. 
The principal argument of the paper was to propose that microsphere buck- 
ling is a dominant contributor to the nonlinear behaviour of the pressure- 
relative volume change curve. 

The computational work carried out in [Shorter et al.l (120081 ) suggests that 
a model incorporating microsphere buckling could successfully predict the 
nonlinearity of the load curve. Therefore here we develop a fundamental 
mathematical model for the loading portion of the pressure-relative volume 
change curve by incorporating volumetric changes and local microsphere shell 
buckling effects. We assume that there is a distribution of shell to radius 
ratio thicknesses and we will suppose that the microspheres are distributed 
dilutely, so that interaction effects between microspheres can be neglected. 
Interaction effects will be considered in future work. 

In order to incorporate the effect of buckling of the microsphere shell, 
we must understand how a spherical shell buckles inside an elastic medium 
under far-field hydrostatic pressure. Although a great deal of classical work 
exists regarding the buckling of spherical shells (where t he imposed pressur e 
is on the sur fa ce of the shell itself ), see for example IWesolowski ( 1967 ): 



Koitei] (Il969h: IWang and Ertepinarl (jl972l ) and more recently IFuI 



( 11998 



Ben Amar and Goriely (|2005[ ). there is a surprising lack of work regarding 



the buckling of shells (of any geometry) that are embedded inside another 
medium. Of specific interest is how the host medium affects the classical 
buckling pressure. 

Initial work into the buckling pressure of a spherical sh ell embedded in an 



unbou nded uniform elast i c med ium has been carried out by lFok and Allwright 



( I2OOII ) and IJones et al.l ( l2008l ). We shall discuss these models and their 



assumptions later on in the paper, particularly that of iFok and Allwright 



( I2OOII ) which is the model that w e shall adopt for buckling here. As de- 
scribed above, [Shorter et al.l (120081 ) also carried out some experimental work 
related to this problem. 

The fundamental objective here then is to determine a model for the 
loading portion of the pressure-relative volume change curve by considering 
a distribution of shell thickness to radius ratio of microspheres which are di- 
lutely dispersed throughout the material. We also introduce nonlinear (finite) 
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elasticity in order to incorporate large deformation of the rubber composite 
in the post-buckling regime. The theory provides a modelling tool to assess 
certain likelihood scenarios. In particular we are able to assess the sensitiv- 
ity of effective properties to changes in specific parameters, e.g. distribution 
of microspheres, nonlinearity and constitutive behaviour of the constituent 
materials that make up the composite and the gas law inside microspheres. 

2. Preliminaries and background 

We consider a composite material with two constituents (or phases as 
we shall term them here) known as the matrix and inclusion phase. The 
matrix phase is a (possibly compressible) homogeneous rubber material and 
the inclusion phase consists of a distribution of thin spherical shells (possibly 
filled with some gas) of initial radius A and shell thickness H. We allow 
for the possibility of a distribution of microsphere shell thickness to radius 
ratios X = H/A. This distribution is governed by a probability distribution 
function F{X). The volume fraction of the inclusion phase is denoted by $. 
We are specifically interested in the problem where the material is subjected 
to an external hydrostatic pressure p. We shall state all pressures relative 
to atmospheric pressure Patm so that upon defining p = p — Patm, p = 
corresponds physically to atmospheric pressure. Similarly, we assume that 
the gas inside the microspheres is also initially at atmospheric pressure so 
that denoting p^^ as the internal hydrostatic pressure we can initially set 
Pin = Pm — Patm = 0. Henceforth all pressures are thus defined relative to 
atmospheric pressure. 

Consider for the moment a single microsphere embedded in an unbounded 
matrix material so that we assume that the pressure is applied in the 'far- 
field'. At a critical far- field pressure Pc, this shell will buckle and the mi- 
crosphere will lose its compressive rigidity for p > p^.. Since we assume that 
we have a distribution of microspheres, each with a different X = H/A, the 
microspheres will buckle successively as the external pressure is continuously 
increased. We illustrate this in Fig. [31 The prediction of the pressure- volume 
relation, given the volume fraction $ of microspheres, the constitutive be- 
haviour of the matrix phase, the elastic properties of the microsphere shell, 
the knowledge of gas internal to the microspheres and the overall distribution 
of X, is clearly a non-trivial problem. The effects of interaction on buckling 
have thus far not been studied and therefore here we consider the case where 
buckling of a microsphere depends only on the far-field hydrostatic pressure p 
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Figure 3: A cube of the microsphere composite material (no scale is implied) is subjected 
to hydrostatic pressure p in the far field with inset figure indicating that the microspheres 
are spherical pre-buckling. This spherical symmetry is retained until the onset of buckling, 
after which the shells deform significantly, losing their rigidity in this state, as indicated in 
the figure inset. Wc do not indicate the precise structure of the shell post-buckling: this 
requires detailed post-buckling analysis. 



and not on the influence of other microspheres. We anticipate therefore that 
this model is vahd for a dilute dispersion of microspheres. We note however, 
that in many homogenization theories, it is often surprising how accurate 



dilute dispersion theories are even in the non-dilute regime (jParnell et al 



2OIOI ). Later work will consider interaction effects in more detail. 



We shall consider each microsphere to have a flxed initial radius A and 
let the shell thickness H vary, so that X = H/A is governed by a probability 
distribution function F{X). Alternatively we could consider H flxed and vary 
A but it transpires that the analysis of the former is more straightforward. 
(Note that for a single microsphere inclusion the scale invariance means that 
varying A ot H for any flxed X must yield the same result.) In some cases 
we need to refer to the middle surface of the shell whose radius we denote by 
A = A — H/2, and the shell thickness to middle radius ratio as X = H/A. 
Also the probability distribution function can therefore be given in terms 
of X, i.e. F{X). With reference to Fig. H] we deflne a (flctitious) radius 
S > A hj the condition $ = (A/S)^ where $ is the prescribed volume 
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Figure 4: A single 'composite sphere' loaded by the hydrostatic far-field pressure p. 



fraction of microspheres. We then consider how the material deforms, and 
how the microsphere buckles, given some hydrostatic pressure p in the far- 
field with this region inside S (containing the microsphere and which we will 
term the composite sphere (CS)) embedded in a purely matrix material. The 
prescription of the radius S allows us to consider the volume change of the 
matrix region under compression. 

We denote by n, fi, E, v the bulk, shear and Young's moduli and Poisson's 
ratio respectively from linear elasticity and we note the relations v = {3k — 
2/i)/(2(3K + fi)) and E = 9K/i/(3K + /x) since later we usually specify fi and 
K. We will make use of the subscripts s and m when we wish to refer to the 
shell and to the matrix medium, respectively. 

Before the microsphere shell buckles (which we shall term the pre-buckling 
stage) we consider the elastic behaviour of both the matrix and microsphere 
shell to be linear. As will be shown, this is reasonable since the shell stiffness 
is significantly higher than that of the matrix phase an d therefore induced 
strains in both media will be small (see also section 7 in iJones et al.l (120081 ) 
for more details). After the microsphere shell buckles (which we term the 
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post-buckling stage) we make the assumption that the shell will lose almost 
all of its rigidity and therefore that the post-buckled microsphere can be 
replaced by a cavity (whilst still ensuring continuity of displacement and 
traction between the matrix and fluid as we shall show later). In this post- 
buckling regime we incorporate nonlinear elastic behaviour by permitting 
large strains and also nonlinear constitutive behaviour. 

In order to justify the linear pre-buckling and nonlinear post-buckling 
assumptions, respectively, we consider the following example. We are here 
interested in understanding how the stiffness of the shell can make the mate- 
rial more rigid as compared with the case when the shell is absent. To this end 



we ca n consider for example the case of a thin glassy shell (see iBaird et al. 



( I1999I )). with surrounding polymeric elastomer composed of a polyurethane 
material. The matrix Young's modu lus can be taken as Em = 3.6 MPa and 



Poisson ratio is typically close to 0.5 (jPiaconu and Dorohoiall2005l ). In terms 



of bulk and shear moduli, we choose the parameter set 

fis = 1.26 GPa, Ks = 2.1 GPa, 

f^m = 1.2 MPa, Km = 4: GPa. (1) 

Note that with this choice, Um = 0.49985 i.e. the matrix is considered essen- 
tially incompressible. Next, let us take a microsphere with shell to radius 
ratio X = 0.01 and for an imposed scaled far-field pressure p/^m, we evalu- 
ate the scaled displacement (with notation reported in section IXT]1 u'^{A)/A 
a.t r = A (the radius at which the displacement is maximum). In Fig. [S]we 
plot this maximum displacement as a function of the imposed pressure when 
the shell is and is not present (left and right on the figure, respectively) and 
in both cases we note that this is predicted by linear elasticity theory. The 
dashed line denotes the critical pressure for this shell to radius ratio, pre - 



dicted by the Fok- AUwright buckling criterion ([S]) f lFok and Allwrightl . l200ll ) . 
When the shell is included, values of u{A)^/A remain small for applied pres- 
sures below this critical value (see the left side of Fig. [5]). On the contrary, 
in the absence of the shell, when reasonable values of the pressure are ap- 
plied, the linear theory is no longer appropriate due to the large values of 
the scaled displacements u^{A)/A in this case (see the right side of Fig. EJ. 
We conclude that in this latter regime we must therefore incorporate full 
nonlinearity in order to permit finite deformations. 

Although we are concerned with matrix materials that ostensibly behave 
incompressibly (typically /im/'^m is small, in particular of the order of 10~^, 
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Figure 5: Plot of the the scaled radial displacement u{A)™/A determined using linear 
elasticity when the shell is (left) and is not (right) included in the model. Associated 
values of X used are X = 0.01 and X = respectively. The dashed line is the critical 
pressure predicted by 



see lOgdenl fll976l )). in section H2] we also calculate the volume change post- 
buckling for a nearly-incompressible theory. As will be shown, it is difficult 
to distinguish the difference between results for the slightly-compressible and 
incompresssible cases, as should be expected. 



3. Pre-buckling behaviour and the buckling model 

In the pre-buckling stage, we consider both shell and matrix phase to 
be compressible linear elastic materials (noting that the matrix is almost 
incompressible) which are perfectly bonded, and we assume that gas resides 
inside the microsphere providing a constant internal pressure p^. We wish 
to determine the total volume change (relative to the initial volume) of the 
material and in order to do this we determine the volume change in the 
composite sphere (CS) when the pressure p is imposed at infinity. 

3.1. Linear elasticity 

Under the assumption of linear isotropic elasticity, the governing equa- 
tions of the corresponding static boundary value problem with no body forces 
are given as follows 

= 0, (2) 
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CTij = XekkSij + 2fieij, (4) 

where aij,eij, and are the components of the stress and strain tensors, 
and the displacements, respectively and we have introduced the Kronecker 
delta tensor 6ij. The matrix material is homogeneous with Lame constants 
A,/i where A = k — 2fi/3. Since the problem is linearly elastic, geometry 
is spherically symmetric and a purely radial stress is applied, then spheri- 
cal symmetry is preserved {ur = Ur{r),ug = = 0). Hence, equation ([2]) 
reduces to a single second order ordinary differential equation which is inde- 
pendent of the Lame moduli. The general solution for the displacement in 
the shell and medium region is therefore of the form 



ulir) 



A,r + §, (5) 



where i = s,m refers to the shell and matrix respectively, and Ai,Bi are 
constants that are fully determined by imposing continuity of displacement 
and radial stress on r = A and the following loading boundary conditions 

<(A-i/) = -p,„ a,"^(r) =-p. (6) 

3.2. Pre-buckling: relative volume change for each CS 

Let us consider a single CS of initial radius 5* and volume V containing 
a microsphere of undeformed radius A and X = H/A. When we increase 
the far-field pressure p {0 < p < Pc), this volume V = (4/3) vrS''^ reduces 
to V = (4/3)7rs^, where s = S + u'^{S) denotes the deformed radius of the 
CS, referring to (E]). The relative volume change, say 6v, occurring in the 
pre-buckling stage is therefore given by 



V \S 

Note that we have assumed the inner pressure inside the microsphere to 
remain constant under loading. This appears to be reasonable since volume 
changes will be small, but will not remain valid in the post-buckling regime 
as we shall consider in section HI 
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3.3. Microsphere buckling 

In this section we discuss the buckhng of a spherical shell inside an un- 
bounded elastic matrix medium, loaded by a fa r-field hydrostatic pressur e 
p. We employ a buckling model introduced by iFok and AUwrightl (200l|). 



Given a distribution of sizes of microspheres inside the material, our aim is 
to determine which of them, for a given imposed pressure p, have buckled 
and which remain unbuckled. 



In iFok and AUwrightl (120011 ) a criterion was derived for the buckling of 
a spherical shell embedded in an elastic material and loaded by a far-field 
hydrostatic pressure under the main assumptions that deformations 
isymmetric and the shell is inextensible. They also neglected the inner gas 
pressure, so in our model we must set pi^ = in the pre-buckling phase. This 
latter simplification is, in fact, not too severe as the displacement is affected 
very little by internal pressure p re-buckling. The a ssumption of axisymmet- 
ric buckling is not a restriction; IWesolowskil ( 119671 ) showed that the critical 
mode number for buckling is the same whether the eigenmode is symmetric 
or not. Also, for glassy shells in the present model, it is easy to show that 
the assumption of inextensibility is consistent with the estimates found for 
the radial and shear stresses. 



Fok and AUwrightl (120011 ) obtained a formula for the critical pressure 



n) in the form 



P{^,n) 21 + lAnf l-Us Em l\f . s^3^ f\V^ I 



n 



(8) 



where X = H/A and n are the shell thickness to middle radius ratio of 
the microsphere and the mode number respectively. Note that in the Fok- 
AUwright approach, n is a natural number greater than 1; dilatational (n=0) 
and rigid-body (n=l) modes are not considered. The functions pi,p2 and ps 
are given by 



p, = [n{n + l)-{l-u,)]/12{l- 
r> -F IF (2n3 - + 3n + 2) 



1/2), p2 = 2/[(n-l)(n + 2)(l + z/,)] 
-i/^(2n3-3n2 + 5n + 2) 



(n - l)2(n + 2) (3n + 2 - 2i/„,(2n + 1)) (1 + i/„) ' 



(9) 



The standard approach would be to specify the shell ratio X and material 
constants h's,h'rn,Es, and substitute these into ([8]) to give the critical buck- 
ling pressure found by minimizing with respect to n (and thus this also gives 
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the corresponding buckling mode n). Here, however we need to approach 
the problem slightly differently since we have a distribution of microsphere 
sizes and we wish to know what the state of that distribution of microspheres 
(buckled/unbuckled) is at a given pressure. It will prove convenient, there- 
fore, to treat n as a parameter as we now explain. Assuming X is for now 
unspecified, by a continuity argument, assuming that n > 1 is a given real 
number we determine the minimum by insisting that dp{X,n)/dn = and 
solve for X. From trivial algebraic considerations it is straightforward to 
show the existence of a minimum for p{X, n) via 

p[{n)X^+p',{n)X + pi,{n) = (10) 

where prime denotes differentiation with respect to argument. The real 
positive root X in f llOp depends on the mode number n, which we de- 
note by Xc(n) where the subscript c refers to critical. Thus, we specify 
n e (l,oo), determine the corresponding Xc from ffTOj) . and then the cor- 
responding pc = p{Xc,n) from (|8]). In this manner we therefore know that 
shells in the range X < Xc are buckled whereas those for X > Xc remain 
unbuckled. 



In Fig. we plot the predictions given by the iFok and AUwrightl (j200l[ ) 
model for the critical pressure Pc as a function of critical size Xc, letting the 
buckling mode parameter n lie in the range [13, 1000], assuming = and 
given the material constants in ([1]). This range of buckling modes gives rise 
to realistic pressures; choosing a lower n corresponds to a higher pressure. 



4. Post-buckling behaviour: nonlinear elastic response 

As we have emphasized before, as we increase the hydrostatic load p grad- 
ually, an increasing number of shells will transition to a buckled state. When 
the shell buckles, there will be some complex modification to the structure 
of the shell and the local matrix medium. Buckling of the shell will re- 
sult in a local loss of rigidity and (initially at least) a macroscopic increase 
in compressibility. Modelling the exact modification to the numerous shell 
structures is a formidable task and hence from a modelling viewpoint, in 
order to model the post-buckling behaviour we shall make the following sim- 
plification. For pressures p > pc, for a CS region, we assume that the spher- 
ical shell region is replaced by a spherical cavity, which at the pressure pc, 
has the same radius as the microsphere at its buckling pressure. The post- 
buckling stage of behaviour will be analyzed under the nonlinear deformation 
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Figure 6: Plot of the critical hydrostatic pressu re Pc/ fJ-m as a function of the critical ratio 
Xc predicted by the Fok and Allwright ( 200l[) model, obtained for n € [13, 1000], when 
we use the material constants specified in ([Ij and we assume pm = 0. 



assumption, because the loss of rigidity of the microsphere will permit finite 
deformations if appreciable pressure p/ fim is applied (note again Fig. |5]). We 
wish to understand how an additional increase in pressure in the CS region, 
past the buckling pressure Pc, decreases the volume of the CS further and 
subsequently we shall derive this influence on the macroscopic volume of the 
composite material. 

From the pre-buckling analysis we can determine exactly the predictions 
of the deformed radii Sc and Oc of the CS (with corresponding initial radii 
S and A respectively) where the subscript c indicates the critical value at 
the buckling pressure Pc- We can therefore begin our nonlinear analysis from 
those values, further increasing the load pressure until we reach the chosen 
load p. This two stage linear-nonlinear approach is however not particu- 
larly appealing; in fact we are able to consider the volume change in the 
post-buckling regime by considering an alternative (nonlinear) problem from 
the outset (i.e. increasing the far-field pressure from zero) that is statically 
equivalent to the linear problem in the pre-buckling regime as we shall now 
show. 

With reference to Fig. [TJ consider a full nonlinear elasticity formulation of 
the deformation associated with an unstressed medium within which resides 
a cavity with the same radius A as the initially undeformed microsphere. We 
consider the deformation of the spherical cavity due to a far-field pressure 
with an additional internal 'shell pressure' denoted by pf^^ which mimics the 
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Increase of external pressure 



Pc P 




Figure 7: To determine the nonlinear deformation post-buckling we use a problem that 
is statically equivalent to the true problem which is drawn schematically here. Above 
is the linear- nonlinear process of deformation and below the fully nonlinear process. The 
presence of the shell is represented by the internal pressure pf^ in the latter, which ensures 
the correct radius when p = pc- 

residual presence of the shell. We choose this pressure by referring to the 
(linear) pre-buckling analysis, considering this to be the maximum pressure 
that the shell exerts on the matrix pre-buckling, i.e. 

= atp = pc, (11) 

such that at the critical pressure we recover to a good approximation the 
linear elastic solution, i.e. to obtain agreement for Oc and Sc- This therefore 
gives the correct starting point for volume change calculations for p > p^. 

We shall consider several assumptions for the nonlinear elastic matrix, and 
we make use of a bar on quantities in the post-buckling regime, in particular 
the undeformed radii A, S and the initial volume V are now denoted by d, 
s, and V respectively. 

We must now determine the deformation of the cavity subject to the im- 
posed external pressure. Since this is spherically symmetric we assume that 
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the deformation is purely radial and we therefore write this radial deforma- 
tion in spherical polar coordinates as 



r{R), 



9 = e, 



:i2) 



where [R, 0, $) are the polar coordinates in the reference configuration and 
(r, 6, (j)) are the polar coordinates in the curre nt configuration , respectively, 
with dr/dR > 0. From the result of Ericksen flErickseru . Il955l ) this deforma- 
tion is not a controllable deformation that is possible in every compressible 
homogeneous and isotropic hyperelastic material. Therefore this inhomo- 
geneous deformation for compressible materials has to be discussed in the 
context of special materials. For example for a special B latz-Ko material an 
analytical solution h as been found in a parametric way (jChung et al.l . Il986 ; 



Horgaru . Il989l . Il995l ). Six other classes of compressible materials have re- 



ceived much attention in the literature where the solution can also be found 



analvticall y (ICarrolj. Il988l . ll991a] Jbl: iMurphyl . Il992l ) . For an overview of such 
results see iHorganI (120011 ) . 

In the incompressible case, thanks to the constraint of incompressibil- 
ity, the radial deformation f|T2|) can be treated in a much more straightfor- 
ward manner and i ndeed is a universal solution (we refer to section 57 of 
Truesdell and Noll! ( 1l992l ) for more details). This deformation satisfies the 
balance equations with zero body force, its equilibrium is supported by suit- 
able surface tractions alone, and is the same for all materials (in the class of 
constrained materials). The choice of strain energy function therefore merely 
dictates the stress field induced by the deformation. 

The special compressible solutions for radial deformation referred to above 
are not well suited to describe the constitutive response of a rubber-like ma- 
terial and therefore we choose not to use these here. Our analysis is therefore 
concerned firstly with purely incompressible materials before we move on to 
describe nearly-incompressible materials in the context of an asymptotic the- 
ory with a small parameter ^u/k <^ 1. In the latter case we use a constitutive 
model proposed in the literature by iHorgan and Murphv (I20091) which was 
based on experimental behaviour (see for example iPenn 

For ease of exposition we have placed details of the theory of nonlinear 



elasticity associated with the subsections to follow in Appendix A 
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4.1. Post-buckling with incompressibility 

The polar components of the deformation gradient associated with ( !T2|) 
are given by 

F = diag(dr/di?, r/R, r/R) (13) 

and for an incompressible material the constraint of incompressibility det F = 
1 states that 

r{R) = {R^ + aY^\ (14) 

where a = — A^. The equilibrium equations in the absence of body forces 
are divT = where T is the Cauchy stress tensor. In the radially symmetric 
case these reduce to the single ordinary differential equation 

^ + - (Trr - Tee) = (15) 
dr r 

where T is derived from a strain energy function (SEF) W, as described in 



Appendix A In this subsection we shall consider two incompressible models, 



the neo-Hookean material with SEF PVnh and Mooney-Rivlin material with 



SEF PFmr; whose forms are given in [Appendix A 



Up to the point of buckling note that we have assumed that the pressure 
inside the microsphere is zero since changes in volume were very small. 
In the post-buckling stage we require the additional pressure pj*^ initially to 
maintain continuity as described in Fig. [7] but we shall also now assume that 
the inner pressure can be non-constant. This is motivated by the fact 
that volume changes can now be large and therefore in this post-buckling 
stage at some point the gas interior to the cavity can be compressed so much 
as to act to stiffen the material. We therefore may anticipate a Boyle's law 
type relation for a massless ideal gas, of the form 

Pin = Pin + Pin (16) 

where p^^ is a constant that accounts for the residual effect of the buckled 
shell and 

p!'n=patm (^i^^y (17) 

In the latter expression the pressure and volume are related through a poly- 
tropic exponent relationship, for a diatomic gas, with exponent rj = 1.4 
where r) is the heat capacity ratio. Note that all pressures are stated relative 
to atmospheric pressure, which motivates the form in (fT7|) . 
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Given a SEF W, it is straightforward to integrate (fT^ and apply the 
traction boundary conditions 



T„.{R) 



R— i-oo 



TrM) 



(18) 



in order to obtain an expression hnking the deformed inner radius a and the 
imposed pressure difference. For the neo-Hookean case we obtain 



P-Pu 



-(- 

2 V a 



whereas for the Mooney-Rivlin model, we find 



A 



5 
2' 



(19) 



P - Pin - Pi, 



- + 7 

2 ' 



-(- 

2 V a 





5 




~ 2 















+ 



• (20) 



1/2 



Note that 7 is a material constant with — l/2<7<l/2, and when 7 
in (120|) . we recover the neo-Hookean solution (fT9l) . 

Using the expressions (fT9l) or (l20l) for a, from (fT4l) we can then obtain 
the (post-buckling) deformed radius of the CS, i.e. s = r{S), predicted by 
the neo-Hookean (1A.3P or Mooney-Rivlin (lA.4p models, respectively. 

4-2. Post-buckling with slight compressibility 

A great deal of work has been presented in the literature regarding the 
constrained theory of elasticity (e.g. incompressible materials) where many 
solutions have been obtained in order to describe approximations to real 
materials. They are approximations because of course no material is in re- 
ality completely incompressible. Numerous constitutive models have been 
proposed in order to model the true behaviour of the material when there 
is a slight deviation from incompressiblity. We assume, as before, that the 
material is homogeneous , isotropi c , and hyperelastic. The early contribution 
to this theory was from ISpencer fll970h and an application was considered 
by iFaulknerl (jl97ll ) who considered the time dependent radial deformation 
of a thick walled spherical shell of almost incompressible material. More 



recent work has been described by Ogden (jOgdenl . 1 19761 . 1 19781 . 119971 ) and 
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HorEfan and Murphy dHorgan and Murphvl . l2007al lbl. l2009h . In order to de- 
scribe the hnearity between pressure and vol ume c h ange (assumed to hold 
for pressures up to 50 MPa) su mmarized by PennI (jl970 ) for natural rub- 
ber, iHorgan and Murphyl (120091 ) derived several different forms of the strain 
energy function and we take the form Wum in flA.7p . 
In the case of slight compressibility, we have 



Krr 



< 1. 



(21) 



We then consider a regular perturbation problem, seeking asymptotic expan- 
sions in powers of e for the relevant solutions, with the results obtained in 
section 14.11 arising as the leading order terms. Thus, we anticipate that the 
leading order deformation will be described by f|T^ and we seek corrections 
to this, associated with the strain energy function (IA.7|1 . We thus derive the 
deformed radii a and s which are slightly modified according to the slight 
compressibility of the matrix. 

We assume that we can write for the deformed radial coordinate 



r{R) = ro{R) + eri(i?) + e^r^iR) + 0(e=^ 



where 



ro(i?) = {R^ + aY^' 
is determined from the incompressible theory (equivalently e — 0). 



(22) 



(23) 



Employing the asymptotic scheme, whose details we provide in Appendix A 
for ease of exposition, we derive the correction term ri{R) in (l22l) as that 
given in ( ]A.16p . This allows us to derive the additional volume change due 
to the compressibility of the matrix medium. 

4-3. Post-buckling: relative volume change for each CS 

For p > pc, let us consider an initial composite sphere of radius S con- 
taining a microsphere of initial size X = H/A. The relative volume change 
for a buckled microsphere is now given by 



Sv 



V -V 
V 



1- ^ 



(24) 



where s = r{S). For an incompressible matrix therefore s is evaluated 
via f lT4|) for Neo-Hookean or Mooney-Rivlin materials. Alternatively, for a 
slightly compressible matrix with strain energy function flA.7p it is evaluated 
via 
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5. Predicted pressure-relative volume change curves for the micro- 
sphere material 

5.1. Total relative volume change for the material 

In sections 13.21 and 14.31 we calculated the relative volume change for each 
composite sphere associated with the pre-buckling and post-buckling stages 
respectively. The choice of the former or the latter depends upon whether the 
value of the applied pressure is below or above the theoretical critical pressure 
required to buckle the microsphere of shell thickness to radius ratio X in 
the composite sphere, according to the Fok-Allwright theory. Our principal 
goal is now to use these two models in order to predict the pressure-relative 
volume change curve for the material as a whole when there is a distribution 
of different shell thicknesses. 

Let us introduce a probability distribution function F{X) which describes 
the distribution of the microsphere shell thickness to radius ratios. Thus we 
impose the far-field hydrostatic pressure p and then use the buckling model in 
section [3l3] to predict the critical Xc below and above which buckling will and 
will not occur respectively. In this way we establish, at each given pressure 
p, the proportion of microspheres that are in a buckled state. Then, in order 
to determine the macroscopic pressure-relative volume change curve, we use 
([7]) and (l24l) for that proportion of microspheres that are in the pre-buckled 
and post-buckled states. 

The relative volume change of the entire material, say (5V, is therefore 
given by the sum of all of the relative volume changes in each composite 
sphere, the distribution of X being accounted for by the probability distri- 
bution function F{X). As such we write 

= I (SviXMX) + 6v{X)x{X)^ F{X) dX. (25) 

where x{^) = 1 ~ x{^) is the indicator function, defined as 

^<-'^{:; 1:13: 

Note that in f l25p we have allowed X to take on all possible values X G [0, 2] 
corresponding to H ^ [0,^]. However, we note that the buckling theory 
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above is applicable to thin shells only and therefore the choice of F is im- 
portant in order for f l25p to give an accurate prediction. In reality the mi- 
crosphere shells are very thin (0(0. 01)) and thus the distribution function F 
must model this, as we now describe. 

5.2. Parameter studies 

Let us now consider how the pressure-relative volume change curve is 
affected by the numerous parameters in the problem and in particular we 
wish to understand the sensitivity of the curve to these parameters. We 
first consider the for m of the probability distribution function F{X). From 



Brazier-Smithl f ll999l ) it appears that the microsphere shell to radius ratio 



distribution can be described well by a Gamma distribution. We thus define 

k x'x'-' 



where > is the shape parameter, Xq > is the mean value (expectation) 
of X and T{k) is the Gamma function evaluated at k. 

Parameter studies will involve choosing values for the elastic properties 
{ks, l^s, ^^m), the initial volume fraction $ of microspheres, the param- 
eters k and Xq in fl27|) and the constitutive model for the nonlinear elastic 
matrix described in section |H Additionally we must decide whether or not to 
incorporate Boyle's law inside the microspheres during compression. This pa- 
rameter set can therefore be chosen in many ways. We select certain cases in 
order to illustrate specific aspects of the model and hence test the sensitivity 
of the results to the particular parameters. 

5.2.1. Influence of nonlinear elastic constitutive model 

We start by fixing the material properties as those given in ([T]) and we 
take an initial volume fraction of microspheres as $ = 0.05. Furthermore, 
within the distribution func tion F{X) giv e n in fl271) we take the shape pa- 



rameter /c = 8 as given in iBrazier-SmithI ( 1l999l ) and we consider a mean 



value Xq = 0.01. Finally, we assume that the gas inside the microspheres 
remains at a constant atmospheric pressure, i.e. pf„ = 0. On the left of 
Fig. [8] we plot the predicted pressure-relative volume change curve for three 
different nonlinear elastic constitutive models in the post-buckling regime: 
Neo-Hookean (dotted), Mooney-Rivlin with 7 = 1/18 (solid) and slightly 
compressible Horgan-Murphy (dashed) as well as a purely linear elastic de- 
formation (dot-dashed). 
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Figure 8: Left: Predicted pressure-relative volume change curve with fixed parameters 
<& = 0.05, Us, Ks, Urn, Km as given in ([T]) and F{X) as in pT]) with = O.Of,fc = 8. 
Nonlinear model for the matrix is taken as neo-Hookean (dotted), Mooney-Rivlin (solid) 
and slightly-compressible Horgan Murphy (dashed) with 7 = 1/18 and for reference we also 
plot the response when the post-buckling regime is determined via linear elasticity (dot- 
dash). Right: Plot of the difference V in the prediction of SV via the alternative nonlinear 
models. V ~ 5Vmb — SVnh (sohd), V = 6Vhm — SVmr (dotted), T> = SVhm ~ 5Vnh 
(dashed) and V = SVlin — SVmr (dot-dashed). 



As can be seen from the left of Fig. |8l all curves are strongly nonlinear, 
exhibiting a softening behaviour under loading after the initial linear elastic 
behaviour pre-buckling at small pressures. As should be expected all four 
curves are initially identical with a slight modification in the post-buckling 
regime where the linear or nonlinear elastic model becomes important. How- 
ever even in this regime the curve is fairly insensitive to the constitutive 
model employed, nonlinear models yielding almost identical results. The lin- 
ear model is fairly close to the nonlinear predictions at these pressures but as 
pressure increases the linear results depart significantly as should be expected 
from this approximate theory. 

In order to better compare the models, on the right of Fig. |8] we plot 
the difference between the relative volume change predicted by the different 
elastic post-buckling models. For example the difference between Horgan- 
Murphy and Mooney-Rivlin is calculated as = 5Vhm — SVmr and anal- 
ogously for the other two possiblities. Given the scale, it is clear that the 
difference between any of the nonlinear models is O(10~^). We reiterate that 
at these pressures the linear elastic model is relatively close to the nonlinear 
models but at higher pressures there is a significant departure. 

Although here we are predominantly interested in p/nm = 0{1) we note 
with reference to Fig. [HI that for both incompressible nonlinear materials, as 
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Figure 9: Prediction of 5V for large values of p/ fim with all parameters as given in Fig. 
|5]and plotted for Neo-Hookean (dotted), Mooney-Rivlin (solid) and slightly-compressible 
Horgan Murphy (dashed) nonlinear models as well as a linear elastic prediction (dot- 
dashed). 



p/yUm oo, SV ^ 0.05, as we should expect in that case (parameters cho- 
sen are those associated with Fig. [8]). The Horgan-Murphy model (dashed) 
predicts a slightly different limit for 6V as shown in Fig. |9] since this takes 
into account the slight compressibility of the matrix. The curve associated 
with the linear elastic model illustrates its unrealistic nature for large defor- 
mations. 

5.2.2. Influence of volume fraction, shell properties and pressure law 

We concluded in section 15.2.11 that the predicted curves are relatively 
insensitive to the nonlinear elastic model employed. As a result of this in- 
sensitivity, let us now model the matrix as an incompressible Mooney-Rivlin 
medium - a standard model for a rubber-like matrix medium. Thus, with 
the solid curve from Fig. [8] as a starting point, also plotted is a reference 
curve in Fig. [TOl and let us vary other parameters in order to assess their 
influence. For each curve we keep all parameters fixed except one control pa- 
rameter in order to assess its particular effect. Thus, the dashed curve in Fig. 
ITU] corresponds to changing only the volume fraction of microspheres from 
$ = 0.05 to $ = 0.1. The dotted curve corresponds to incorporating Boyle's 
law ( fT7|) for the gas interior to the microsphere instead of constant pressure. 
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The thick dashed curve corresponds to softer shell material properties, i.e. 

/i, = 0.126 GPa, = 0.21GPa, (28) 

and the dot-dashed curve is associated with slightly more compressible matrix 
material properties (for the linear elastic pre-buckling portion of the curve), 
i.e. 

= 1.2 MPa, Km = 0.4 GPa. (29) 

These correspond to a Poisson's ratio of Um = 0.4985. 

Let us assess each one of these in turn. Increasing the microsphere volume 
fraction, whilst keeping their distribution fixed has the expected effect: the 
curve remains qualitatively similar, becoming relatively softer in the nonlin- 
ear regime. Incorporating Boyle's law should add some stiffness in the post- 
buckling regime and this can be seen in the figure; but its effect is rather 
modest. The dotted and solid curves are identical until the post-buckling 
effects become important, around p/fim = 0.2, and then increased pressure 
interior to the microsphere does yield a small additional stiffness. Softer 
shell properties are expected to have a larger effect in the transition region 
from pre to post buckling as the shells will clearly buckle at lower pressures. 
This can clearly be seen in the figure; an order of magnitude change to the 
properties has modified the curve significantly. However, the insensitivity 
post-buckling can be seen by virtue of this curve and the solid curve remain- 
ing parallel in this regime. Finally, the small decrease in matrix Poisson's 
ratio corresponding to the dot-dashed curve yields the expected effect; the 
material becomes slightly softer in the linear region, recovering an identical 
nonlinear response to the reference Mooney-Rivlin case in the post-buckling 
region. 

As perhaps should be expected, there is a great sensitivity to the mate- 
rial properties of the shell, but the model is relatively insensitive to other 
parameters. 

5.2.3. Influence of probability distribution function parameters 

Let us now take the Mooney-Rivlin reference curve (solid) as plotted in 
Fig. [8] and vary the distribution function parameters in F from those of 
the reference material Xq = 0.01 and k = 8. In Fig. [11] we plot the distribu- 
tion function (left) and corresponding pressure-relative volume change curves 
(right) whilst keeping k = 8 fixed and varying Xq (top) and then keeping 
Xq = 0.01 fixed and varying k (bottom). As perhaps should be expected 
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Figure 10: Parameter study for the predicted pressure-relative volume change curve. The 
sohd curve corresponds to the parameters assumed in Fig. [8] with a Mooney-Rivhn matrix 
medium. We then vary one specific aspect of the model to assess sensitivity. Different 
curves correspond to: increasing microsphere volume fraction from $ = 0.05 to $ = 0.1 
(dashed); incorporating Boyle's law P7)) for the gas interior to the microsphere instead 
of constant pressure (dotted); softer shell properties given by (|28p (thick dashed); slightly 
more compressible matrix phase given by (|29p (dot-dashed). 

there is great sensitivity to the average shell thickness to radius ratio. This 
can be seen by the large modifications to the pressure-volume curves in the 
top-right figure. The main effect of varying k is the later (in terms of higher 
pressure) softening of the material, although its influence is less marked than 
variation in Xq. This information is useful from the viewpoint of knowing 
the correct type and distribution of microspheres to use in the composite. 

Finally, we illustrate how the 'kink' in the load curve is associated pri- 
marily with the distribution function of the shell thicknesses. Let us choose 
Xq = 0.01 and successively increase k which has the effect of tending the 
distribution function towards a dirac delta function as can be seen in Fig. 
[T2] where we take k = 8 (solid), k = 50 (dotted) and the limit as /c — )■ oo 
(dashed). The limiting case when there is only one size of microsphere shell in 
the composite {k — )■ oo) is reflected in the shape of the pressure- volume curve 
illustrated on the right of Fig. [T21 manifested by a discontinuity in deriva- 
tive of the curve. When k is made flnite and reduced the curve becomes 
progressively smoother. 
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Figure 11: Plots of the distribution funetion (left) against shell thickness to radius ratio 
and corresponding pressure-relative volume change curves (right). The parameter fc = 8 
is held fixed and we vary Xq (top) whereas wc keep Xq = 0.01 fixed and vary k (bottom). 
Specifically we take Xq = 0.01, 0.02 and 0.005 (solid, dotted and dashed respectively) (top) 
and k = 8, 15 and 30 (solid, dotted and dashed respectively) (bottom). 
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Figure 12: Illustrating how the shape of the load curve is modified in the limit A: — > oo. We 
choose Xq = 0.01 and take fc = 8 (sohd). A: = 50 (dotted) and the limit as /c — > oo (dashed). 
The limiting case when there is only one size of microsphere shell in the composite (fc — )■ oo) 
is reflected in the shape of the pressure-volume curve illustrated on the right, i.e. the 
appearance of a discontinuity in derivative of the curve. The curve becomes smoother in 
this region as fc is progressively reduced. 
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6. Conclusions 



We have presented a model that predicts the nonlinear pressure-relative 
volume change loading curve associated with a microsphere elastomeric com- 
posite material. The nonlinearity is induced by several mechanisms: (i) 
incorporating a distribution of sizes of spherical shells which thus buckle 
successively according to the applied load, (ii) modelling the post-buckling 
behaviour of the matrix as a nonlinear elastic material, (iii) incorporating 
Boyle's law for the pressure interior to the microsphere in the post-buckling 
regime. 

In this initial study we have neglected any interaction between micro- 
spheres both in terms of the buckling analysis and the determination of the 
change in volume of the composite. Therefore we anticipate that the model 
is valid only for low volume fractions of microspheres. 

Parameter studies reveal that, although it appears important to include 
nonlinear behaviour in the post-buckling stage, the curves are largely insen- 
sitive to the chosen nonlinear elastic model. However the curves are particu- 
larly sensitive to the properties of the microsphere, including shell properties 
and the distribution of shell thickness, particularly the choice of mean shell 
thickness to radius ratio Xq. Furthermore, the shape of the 'kink' in the 
load curve is associated primarily with the distribution function of the shell 
thicknesses. The smaller k is, the smoother the transition to the nonlinear 
post-buckling state. 

We have also seen that when Boyle's law is incorporated in the post- 
buckling regime (the influence of a gas inside the shell pre-buckling is negli- 
gible as so is omitted in the model) there is competition between the softening 
of the material due to microsphere buckling and stiffening due to Boyle's law. 
It transpired that the latter contribution is small, as is visible in the tran- 
sition from pre to post buckling in Fig. [TOl The near-incompressible nature 
of the matrix medium dictates that eventually the response of the composite 
will change, as the pressure increases, from an initial softening material to 
one which stiffens for large pressures {p/f^m > 0{1)) as is seen in Fig. O 

Note that in order to solve the post-buckling problem of a single spheri- 
cal cavity embedded in a unbounded medium subject to inner and external 
hydrostatic pressure we used the theory of almost-incompressible materials, 
posing an asymptoti c expa nsion for the Horgan- Murphy model (1A.7P (see 
Horgan and Murphvl tOO^ ). 



Follow-on work will consider the accuracy of the buckling model, by com- 
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paring with alternative (e.g. iJones et al. and new models. We shall 



also consider the effect of the interaction of microspheres on buckling. This 
is clearly important when volume fractions of microspheres become larger, a 
common case in practice. 
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Appendix A. Nonlinear elasticity theory 

Appendix A.l. Incompressible theory 

Given the deformation gradient F, we compute the physical components 
of the left Cauchy-Green strain tensor B = FF'^ from f|T3|) and find its first 
three principal invariants from 

h = tiB, h = ^[(trB)2 - tiB^], h = detB. (A.l) 

For an isotropic incompressible hyperelastic solid, the Cauchy stress tensor 
T is then related to the strain via 

T=-qI + 2WiB - 2W2B~\ (A.2) 

where q is the Lagrange multiplier introduced by the incompressibility con- 
straint, W{Ii,l2) is the strain-energy density and Wi = dW/dli. 

In the incompressible case, in order to model the nonlinear constitutive 
response of the matrix material, let us consider two alternative strain energy 
functions, the so-called neo-Hookean and Mooney-Rivlin models: 

W^NH = ^(/i-3), (A.3) 
Wmk = \{\ + i) I'mih - 3) + ^ - 7) f^mih - 3), (A.4) 
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where fim > is the shear modulus of the matrix (as introduced above) for 
infinitesimal deformations and 7 is a non-dimensional constant in the range 
-1/2 < 7 < 1/2. 

Appendix A. 2. Slightly- compressible theory 

In the case of generally compressible materials, the Cauchy stress can be 
written as 

T = /3o/ + AB + /3_iB-\ (A.5) 



where 



dW 



dh 



dW 



2 dW 



(A.6) 



-2 a//. 



dW 
'dh' 



A variety of strai n energy functions can be p roposed, but one in particular 
was described by iHorgan and Murphyl (120091 ) in the form 



HM 



^ U + 7 ) (/i - ^ll") + ^ ( ^ - 7 ) (/2 - 3/3^^) + '^{l 



.2/3x 



K. 



■m , jl/2 



(A.7) 

where 7 is an arbitrary constant, is the infinitesimal shear modulus and 
t^rn is the infinitesimal bulk modulus. The latter two material parameters are 
usually related for an almost incompressible (slightly compressible) material 
through the additional assumption 



< 1. 



(A.8) 



Note from ( 1A.7P that Whm can be considered as a slightly compressible gen- 
eralization of the incompressible Mooney-Rivlin strain energy function (1A.4P 
since 

W^hm(/i, /2, /3 = 1) = W^mr(/i, h). (A.9) 

We proceed to solve the boundary value problem by posing the radial 
displacement as a regular asymptotic expansion in powers of e, i.e. 



r{R) = ro{R) + eri{R) + e^r2{R) + 0{e^), 



(A.IO) 
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where 

ro{R) = {R' + af^ (A.ll) 

is determined from the incompressible theory (equivalently e — )■ 0). The 
correction terms ri{R) and r2{R) are the same order of magnitude as ro(-R) — 
R. Making use of the unconstrained theory our goal is to determine ri{R) 
and in order to do this it is necessary to retain terms of O(e^) in (lA.lOp . 

Let us expand the Cauchy stress tensor T, whose form for the compress- 
ible problem is given in (lA.Sp . in the form 

T = To + eTi + 0(e2). (A.12) 



Since f3o in (IA.6P invol ves the t e rm h dW/dl^, the hydrostatic part of Tq 
involves the term ri (see ISpencerl ( 19701 )) which is then determined in a man- 



ner described below. We write all other equations in the form of asymptotic 
expansions. In particular we note that 

divT = divTo + e divTi + 0{e^), (A. 13) 

and the traction boundary conditions must be written as 

Trr{R)\R^oo = -p + 0e + O{e^), Trr{A) = -p',^-p',^ + 0e + O{e^), (A.14) 

noting that they give inhomogeneous conditions at leading order and homo- 
geneous conditions at higher orders. Using these we can then equate terms 
at each order in e and solve the resulting problems. 

Using the boundary condition (lA.14P i to leading order and the equilib- 
rium equation divTo = we obtain the general expression for To^r, 

TUR) = - 2 r ^'fp)T^rr{p)-Toeeip) ^^^^^^ 
Jr dp ro(p) 



Imposing the leading order boundary condition (lA.Mp ? then allows us to 



recover the relationship (120|) . i.e. the result associated with the incompressible 
solution (fT4|) . 

Next equating (IA.15I) with the zeroth order term Torr found in (1A.12I1 . 
allows us to find a linear first order ordinary differential equation for ri, 
whose solution is given explicitly by 

2i?V2(-l + 27) + R\l + 27) + R^ro f 1 - 67 - ^) + 4Ciro 

ri{R) = ^ '-^ , 

4ro 

(A.16) 
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where Ci is an integration constant to be determined. In order to determine 
Ci we must continue our analysis to the next order i.e. we solve divTi = 0. 
The first order balance equation is a linear second order differential equation 
for the unknown function r2{R) whose solution we do not reproduce here (for 
brevity, but also because we are primarily interested in ri (/?)). The function 
r2{R) involves two constants of integration which we shall call C2 and C3. 
It is straightforward to determine that only one of these constants, say C2, 
is involved in the component Ti^r of the stress, in addition to Ci. Using 
the 0(e) boundary conditions (1A.14P for the stress Ti, it is then possible to 
determine both Ci and C2. In particular since we are interested in only the 
first order correction to Tq we state only Ci which takes the form 



where 



480a (A3 + a3) {a^{l - 27) + ^2(1 + 27)) /x. 



(A. 17) 



A(40p(7 
240^2 d 



67) - 3(83 
27)Vm, 



47(437 -57))/iJ, 



-20^3(27 - l)(4p + 3(67 
120A^ (472 - 1) /i^, 
-10A^{p{A + 87) + 3(5 + 47(97 



0^4 

«2 

ai = -40A^(1 + 27)(4p + 3(67 
ao = 135^9(1 + 27) Vm- 
In contrast to iFaulknei] (Il97lh . here the const ant C^ depends on the bound- 



16A^(10p(27 - 1) 
mA' (472 - 1) /i, 
-40A^ 

19/ 



9/im + 67(167 - 9)/i^ 



l)/^m), 



ary conditions via the presence of p and a. iFaulknerl ( Il97ll ) does not im- 
pose boundary conditions in the form flA.14p working instead with rather 
non-physical boundary conditions for the slightly compressible part of the 
deformation. 
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